Effects of pressure on the superconducting properties of magnesium diboride 
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We discuss the effects of hydrostatic pressure on the superconducting properties of MgB2 within 
the framework of Eliashberg theory. By considering the pressure dependences of all parameters 
appearing in the McMillan formula, we show that the calculated pressure derivative of T c as well as 
the variation of T c with pressure are in good agreement with recent measurements. The pressure 
dependences of the energy gap Ao, the effective interaction strength N(Ef)v, the critical magnetic 
field H c (0), and the electronic specific heat coefficient 7 are also predicted for this system. A 
comparison of pressure effect in non-transition elements clearly suggests that MgB2 is an electron- 
phonon mediated superconductor. 

PACS numbers: 74.62.Fj, 74.70.Ad 



I. INTRODUCTION 

The recent discovery of superconductivity in MgB2 
(Ref. |l|) has attracted considerable interest in the study 
of this material, both to understand the mechanism of su- 
perconductivity and to explore other properties of MgB2 
and related materials. The high transition temperature 
T c w 40 K in this material offers other possibility for 
finding high-T c superconductivity in some binary inter- 
metallic compounds besides cuprates and Ceo-based com- 
pounds. MeanwJhile, the high critical currents observed in 
MgB 2 thin filmsB and wiresa reveal that MgB 2 belongs to 
a new class of low-cost, high-performance superconduct- 
ing materials for magnets and electronic applications. 

Measurements of the isotope effect and of the influ- 
ence of pressure on the transition temperature and criti- 
cal field of superconductors yield information on the in- 
teraction causing superconductivity. Indeed, the pressure 
(or volume) and the mass number would seem to be the 
only variables whose effect might be capable of imme- 
diate theoretical interpretation. By observing how pres- 
sure changes the parameters of the lattice in the normal 
state and in the superconducting state, and comparing 
the measurements with the theoretical predictions, one 
can. test the validity of some theoretical models. Olsen et 
aln have shown that the volume (V) dependence of ef- 
fective interaction N(Ep)v, din. N (E p)v / dh\V , can be 
scaled well with the deviation £ from the full isotope ef- 
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fect where £ is defined by T c m. M' 
ducting metals. Bud'ko et aln and Hinks et aln reported 
a sizeable isotope effect for B 26(3) or 0.30(1)) in 

newly discovered superconductor MgB 2 . Although the 
total isotope coefficient a=0.32(l) (Ref. ||) is smaller than 
the canonical BCS value of 0.5, it is the same as that in 
Cd (Ref. 0). The isotope effect along withxiLher measur^ 
ments such as inelastic neutron sca tteripg ,mil tunneling,E£l 
NMR (Ref. 0), and specific heatBBEil confirmed that 
MgB 2 is an electron-phonon mediated s-wave supercon- 
ductor. 

Soon after the discovery of superconductivity in Mb 
the effect of pressure on T c was studied by two groupsc 
by resistivity or ac susceptibility measurements. Both 



groups observed a decrease of T c with increasing pressure, 
with initial pressure derivative dT c /dP of -0.8 K/GPa 
(Ref. HI) or -1.6 K/GRa_(Ref. |TJ), respectively. More- 
over, Monteverde et al£3 found that the superconduc- 
tivity is not destroyed by applying high pressure up to 
25 GPa, at which point T c is as high as 21 K. A some- 
what larger diU/dP of -2.0 K/GPa was late reported 
by Saito et alS-H from high-pressure resistivity measure- 
ments. Using a He-gas apparatus, Tomita et al£3 de- 
termined a dT c /dP of -1.11 K/GPa under pure hydro- 
static pressure conditions. In order to find the reason 
why the reported values of dTJ-rtP are different among 
different groups, Lorenz et al.c3 carried out high pres- 
sure experiments on MgB 2 samples with different T c 's 
at ambient pressure and different pressure media. T c 
was found to decrease linearly over the whole pressure 
range (0-1 GPa). In the He environment, the two sam- 
ples with the initial T c =39.2 and 37.4 K yield the pres- 
sure derivatives of -1.07 and -1.45 K/GPa, respectively. 
Tip, former is obviously very close to that of Tptnita et 
aZ.ll3 The latter approaches their previous datajla which 
was obtained by using the Fluorinert FC77 as pressure 
medium. They therefore concluded that the variation in 
the value of dT c /dP by various groups results from the 
differences in sample preparation conditions. The value 
of dT c /dP K/GPa is then confirmed to give the 

true hydrostatic pressure dependence of T c in MgB 2 . 

Two theoretical models have been tried to describe the 
systematics of the behavior of T c under pressure in MgEL_ 
Based on the theory of hole superconductivity, HirschEj 
predicted an increased T c with the decrease of B-B in- 
traplane distance under the application of in-plane bi- 
axial pressure. However, this prediction has not been 
confirmed experimentally yet. No uniaxial pressure mea- 
surement was reported due to the extreme difficulty in 
growing MgB 2 single crystal. The experiments of hydro- 
static pressure effect on T c do not particularly support 
this theory provided that no charge transfer between the 
Mg and B layers occurs. A4te«w4ittely, the experimental 
results have-been analyzedi3li§ElEl by using the McMil- 
lan formulae^ derived from Eliashberg theory,Ej support- 
ing electron-phonon mediated superconductivity. Inter- 
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estingly, Loa and SyasseiH analyzed the pressure ef- 
fect on T c from McMillan formula on the basis of their 
calculated elastic and electronic structure data. As- 
suming that the electron-ion matrix element / is pres- 
sure independent, they found that the pressure effect 
on T c is in good agreement with experimental data by 
using a lattice Gruneisen parameter jq=1. These as- 
sumptions, however, deserve some refinements. Recent 
band structure calculations suggeajLihat MgB 2 is a tra- 
ditional sp metal superconductor .EjoLj The pressure de- 
pendence of I has long been an interesting i ssue , of t he. 
research of pressure effect in simple sp metals ,r 8 tr HP'ti 3 4c 2 ] 
Ziman's calculation of the electron-phonon interaction 
leads to < IjL_>oc N(Ef)~ 2 at least in the limit of long 
wavelengths. E3 This then indicated that the considera- 
tion of the pressure dependence of / would be important 
for better understanding the superconducting properties 
of MgEL under pressure. On the other hand, it has been 
foundEj that the choice of lattice Gruneisen parameter 
7G is crucially important in explaining both the magni- 
tude and the sign of the pressure derivative of T c when 
using McMillan formula. The value of 7g=1 m the cal- 
culation of Loa and..Syassen is obviously lower than those 
reported recently.oE3 The pressure dependence of the ef- 
fective electron-electron Coulomb repulsion appearing 
in the McMillan formula is usually neglected in previous 
studies due to the assumption of the small change of \i* 
compared with that of the electron-phonon coupling pa- 
rameter A (Ref. |3^). However, the magnitude of fi* is 
also of interest in connection with the possibility.. that 
superconductivity may be destroyed by pressure £30 It 
was argued that the pressure dependence of /i* makes a 
significant contribution to the behavior of T c uxLdep-yery 
high pressures and must be handled carefullyE3't3a 

In this paper we discuss the pressure dependences 
of some interested superconducting properties in MgLV 
The outline of this paper is as follows: In Sec. II we 
presented a theoretical approach for pressure effects on 
the superconducting properties in the simple sp metals 
superconductors. Section III contains the theoretical re- 
sults obtained and a comparison with experiments for 
MgB2. We draw conclusions in Sec. IV. 



II. THEORETICAL FORMULATION 

For our purposes, the relation between T c and micro- 
scopic pas^meters is given adequately by the McMillan 
equation^ 



O, 



1.45 



exp 



1.04(1 + A) 
A -//*(! + 0.62 A) 



(1) 



which relates T c to the electron-phonon coupling param- 
eter A, the Coulomb repulsion strength (i* , and a tem- 
perature Qd characteristic of the phonons. 

Considering the variations of Qp>, A, and /z* with 
pressure or volume and introducing parameters if = 



d\nX/d\nV and <f> = <91n/i* /dhiV : we can get the pres- 
sure coefficient of T r 



dlnT c _ 7G I.04A(1 + 0.38^*) <p 
dP ~ ~B~ ~ [A-a**(1 + 0.62A)] 2 B^ 
f.04/i*(I + A)(f + 0.62A) <j) 
[A-^*(l + 0.62A)p ' 



+ 



(2) 



where B = l/nv = —dP/dlnV is the bulk modulus 
and Qd is assumed to be proportional to < uj 2 > 1 / 2 
and 7q = —d In < uj 2 > 1 / 2 /dlnV being the effective 
Gruneisen parameter. 

It is well known that the usual BCS result for the en- 
ergy gap can be expressed byc3 



A = 29d exp 



N{E F )v 



(3) 



Where N(Ep) is the electronic density of states at the 
Fermi energy Ep and v is the pairing potential arising 
from the electron-phonon interaction. If we renormalizc 
the Morel- Anderson resultcil by introducing the renor- 
malization parameter Z n (0) = 1 + A into their analysis, 
the effective interaction strength N(Ep)v can be rewrit- 
ten asca 



N{E F )v 



X 



I + A 



(4) 



The logarithmic volume derivative of N{Ep)v is then 
given by 



d\nN{E F )v 



A(l + 



-w- 



(5) 



d\nV {X- n*){l + X)~ r A — a* 

Considering the experimental observations of the pres=. 
sure dependence of the energy gap of superconductor,CJ 
we differentiate Eq. (3) with respect to pressure 



dlnAo tg 
B 



dP 



I 



A(1+m*) M*(l+A) 



(A-/i*) 2 "^ (A-/i*) 2 



(6) 



BCS expression for the critical field H c at absolute zero 
temperature isc3 



H c (0) 2 
8tt 



2N{E F )tu 2 exp 



N(E F )v 



(7) 



Differentiating Eq. (7) with respect to the pressure, one 
obtains an expression of the pressure coefficient of H C (Q), 



dlnff c (0) dlnA 



dP 



dP 



2^ + ^ ' ,S1 



f 

2B 



where 7^ = dlnN(Ep)/dln V. 

The expressions for 7c ip, and <j> can be integrated to 
give 



e D (v) = e D (o) [v/Vor^ 

X(V) = X(0)[V/V f 

f(y) = fMlv/Vo]* . 



(9) 
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Here V and Vq are the unit cell volumes under the 
applied pressure and at ambient pressure, respectively. 
These two volumes can be related according to the first- 
order Murnaghan equation of state V(P) — U(0)(1 + 
B' P/B )~ 1/B '°. The Eq. (9) is then rewritten as 



HP) 
M *(p) 



e D (o) 

A(0) 
M *(0) 



B' P 



1G/B' 



1 



B 

B' P~\ ~ v,B ' a 



(10) 
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B 

B'Pl " i,IB ' a 



Bo 



From Eqs. (1) and (10) we arrive at the expression for 
the pressure dependence of T c 



T c (P)=T c [e D (P),X(P)^*(P)} 



(11) 



Knowing Bo, B' , jq, 7^, 0, and tp, one can evaluate 
the pressure effects on the superconducting properties, 
especially the behavior of T c under pressure. Bo and 
B' Q can be obtained from the compressibility data deter- 
mined by neutron or synchrotron x-ray diffractions. A 
direct experimental determination-af 70 can be made by 
measuring-electron tunnelingE3EJc3 or inelastic neutron 
scatteringeJ under high pressure. In general, for met- 
als in which different techniques yield similar Griineisen 
constants, a good approximation to 7^ is provided by the 
room temperature value determined from the Griineisen 
equation 



1G 



KyCp 



where ay is the volume coefficient of thermal expansion, 
V m is the molar volume, and C p is the molar heat capac- 
ity at constant pressure. The approximation for 7^ of 
Slater is .derived from the pressure derivative of the bulk 
modulusa 



g = ^o _l_ _2_ 1 Vd 2 P/dV 2 
lG = 2 _ 6 ~ ~3 ~ 2 



(13) 



dP/dV ' 

The formula for /i* due to Morel and Andcrson0 used 
here is 

M* = 1 , * 7- , ■ (14) 



1 + n\u(E F /uj ph ) 

with (i = 0.51n[(l + a 2 ) /a 2 } and a 2 = Tre 2 N(E F )/k 2 F , 
from which we evaluate the volume dependence of fi* as 



1 



3~ 7G 



2fi 2 



N 



(15) 



Here the variation of k F with volume has been calculated 
from the fundamental definition kp = (S^Z/V) 1 / 3 with 
Z the valency. Unfortunately, to the best of our knowl- 
edge 7jv has never measured directly for any supercon- 
ductor in the case of free electron gas it would have a 



value of 2/3. Using the expression given by MigdalS for 
the electronic specific heat coefficient 7 one obtains for 
the electronic Griineisen parameter 



d In 7 
dlnV 



A 



In 



1 + A 



(16) 



The electronic Griineisen parameter 7 e is usually-deduced 
from measurements through the simple relatiorJ£3 



7e 



a e Vm 

KyC e 



(17) 



Here a e is the contribution to the expansion coefficient 
from the electrons at lower temperatures and C e is the 
electronic heat capacity. A theoretical estimate of j e 
can also be given from the measurement of the volume 
dependence of the orbitaLsusceptibilityEiJ or from band- 
structure considerationsEj 

The electron-phonon coupling parameter A can be ex- 
pressed as 



A = 



NjEp) < I 2 > = 
M < uj 2 > 



V 



M <oj 2 > 



(18) 



where < I 2 > is the mean-square electron-ion matrix 
element and M the ionic mass. The McMillan-Hopfield 
parameter 77 (or N(Ep) < I 2 >) has been regarded as a 
local "chemical" property of an atom in a crystal. Allen 
and Dynesu pointed out that r\ is the most significant 
single parameter in understanding the origin of the high 
T c of conventional superconductors. For strong coupling 
systems, variation in rj is more important than variation 
(12) of < lo 2 > in causing T c to change. Softening of < u 2 > 
often does enhance T c , but very high T c should be caused 
more by large r\ than by small < uj 2 >. 

The logarithmic volume derivative of A, <p, is then ob- 
tained 



ip = 



dhiri 
<91nU 



27g = 5 + 2 7g . 



(19) 



In order to understand how the electronic contribution 
Tj = N(Ep) < I 2 > pu'ies with volume, we use the 
Gaspari-Gyoriffy theorja for rj, i.e., 



kp 



TT 2 N(Ep) 



E 



2(Z + l)sin 2 (^ +1 -Si)NiN t+1 



( 2 °) 

where Ni is the Ith angular momentum component of the 
density of states, is the Ith component of the single 
scatterer density of states evaluated at Ep, and Si the 
phase shift. 

For simple metals, the scatterers are assumed to be 
weak. We can take Ni — and approximate sin 2 (<5; + i — 
5i) by (61+1 — Si) 2 . Eq. (20) is rewritten as 



kp 



TT 2 N(Ep) 



£2(1 + 1)$ 



1+1 



^) 2 



(21) 
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This expression iSpidentical to the pseudopotential for- 
mula of McMillanHj'LJ Assuming that the phase shift Si 
does not vary very much under pressure for simple sp 
superconductors, we then obtain 

S = -lN-l ■ (22) 

The form in Eq. (22) is the same as that of Baryakhtar 
and Makarov,E3 who used the constant of-the electron- 
phonon interaction of Frohlich and Mitra£3 The expres- 
sion is an improvement over the expressions of 5 , =j9-iand 
S = —4/3 obtained by Olsen et alS3 and Seiden,c3 re- 
spectively. It is interesting to notice that substitution 
of 7A r=2/3 into Eq. (22) yields iL= -4/3. Eq. (22) re- 
duces to the expression of Seiden,E3 who modified McMil- 
lan's expression for A somewhat by considering the ef- 
fects of a real lattice spectrum as opposed to the Jel- 
lium model. Since the electronic Griineisen parameter j e 
usually varies among diffc ropt. metals even in the simple 
non-transition elements, Qe3e3 we believe that Eq. (22) 
should provide a more reasonable value of S compared 
with Seiden's formula. 



0.0 




_ 2 i , , , , , , , 1 

1.0 1.5 2.0 2.5 3.0 

Yg 



FIG. 1: Pressure derivative of T c as a function of the lattice 
Griineisen parameter 7g in MgB2. The circles show the cal- 
culation from the four different values of 7 g=1.2, 1.83, 2.3, 
and 2.9. 



III. RESULTS AND DISCUSSION 

Using the experimental value of T c =39.25 K (Refs. 
|l9| ), and the theoretical estimates of A=0.87 and /Lt*=0.10 
(Ref. |f), we got 9d=860 K from Eq. (1) for MgB 2 . 
We believe that all these parameters, which will enter 
our calculations, are reliable. Foe. example, the inelastic 
neutron scattering measurementsa provide an estimate of 
A ~0.9, which is close to that we used. The calculated 
value of 0£>=86O K is in the range from 746 to IQ. 
determined from the specific heat measurements.! 12 !! 13 !!!. 

We took the structural parameters I?o=147.2 GPa and 
B' =4 from the measurements under the pure hydrostatic 
pressures up to 0.62 GPa (Ref. |59] ) and under high pres- 
sures up to 15 GPa (Refs. 22 33), respectively. To our 
knowledge no inelastic neutron scattering or tunneling 
data exist for MgB 2 under hydrostatic pressure. Wc 
have to use Eq. (12) or (13) for estimating the lattice 
Griineisen. parameter jq. The measurements of heat 
capacity^ give a C p of 47.80 J/(K mol) at T=298.16 
K. U m = 1.75xl0~ 5 m 3 /mol, k v =6.79x 10~ 12 Pa" 1 , and 
ay=2.22xl0 -5 p-K" 1 can be drawn from the neutron 
diffraction data.H We therefore obtained 7g=1.2 by us- 
ing Eq. (12). Based on the first-order Murnaghan equa- 
tion for V(P) and the Slater expression of Eq. (13), 
we got a somewhat larger jg of 1.83 compared to that 
from Eq. (12). For most simple metals, there is no 
much difference between the room temperature lattice 
Griineisen parameter given through Eq. (12) and the 
Slater relation.E3E2l It was found that the Slater expres- 
sion usually pGan yield the reasonable values of for 
most metalsJ^ 2 ! The only uncertainty entering Eq. (12) 
in our calculation comes from the indirect measurements 
of the linear coefficients of thermal expansion.LJ Roundy 



et alE reported a value of 7g ~2.3 from ab initio calcu- 
lations, which is close to our calculated fy£ according to 
Eq. (13). Meanwhile, Goncharov et alE5 determined a 
large E^ g mode Griineisen parameter of 2.9±0.3 from the 
measurements of Raman spectra under pressure. This 
value is obviously larger than those derived from Eqs. 
(12) and (13). 

In the calculations of electronic density of states, Loa 
and SyassanEil found that N(Ep) decreases with pres- 
sure at a rate of d\nN(E F )/dP =-3.1xl0~ 3 GPa" 1 . 
Combining this calculated value and experimental value 
of Bo, we got 7at=0.46. The volume dependence of fj,* 
is then derived from (15) once having the values of 7g 
and /i. For simple sp metals, a 2 has a typical value of 
0.4 (Ref. |4l]), which yields fj, =0.63. The volume depen- 
dence of A, ip, is therefore determined from Eqs. (19) and 
(22). Using 7at=0.46, we have £=-1.13 for MgB 2 , which 
is smaller in magnitude than -4/3 in Seiden's formula, 
for simple metalsEj For the transition metals, HopfieldE 3 ] 
commented that S is a relatively constant quantity with a 
value of about -3.5. The values of «S=-3.5~-3.1 obtained 
by inverting the measured dT c /dP for YNi2- a ;M. I ;B2C 
(M=Co and Cu) (Ref. [34j) are comparable to that of the 
transition metals, but are larger in magnitude compared 
to that of MgB 2 . 

With the parameters determined above, we have cal- 
culated the pressure derivatives of T c for MgB 2 by using 
Eq. (2). In Fig. 1 we plotted dT c /dP as a function of 
7g in the interested range. It is interesting to note that 
7G plays a predominant role for the pressure effect of 
T c . For the four different 7g's considered here dT c /dP 
are negative. The values obtained from 7g=1.83 and 
2.3 are -0.78 and -1.12 K/GPa, respectively. These are 
close to the hydrostatic pressure value of -1.1 K/GPa 
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FIG. 2: Variation of T c with pressure in the region of to 1.0 
GPa of MgB2 for 7g = 1.83 and 2.3, respectively. The circles 
and squares represent the hydrostatic presspje experimental 
datetntaken from the works of Tomita et al&3 and Lorenz et 
al.f3 respectively. 



(Refs. |T^ , |l9|) . Thus the hydrostatic pressure results can 
be reproduced in terms of our present model by using 
the values of 7g obtained from cither the Slater relation 
or ab initio calculation. It is difficult to obtain the mea- 
sured results by using 7g=1 as suggested by Loa and 
SyassenrJ We noticed that a jg of 2.27 is necessary so 
as to account for the pressure effect on T c for MgB2. As 
emphasized above, all quantities entering Eq. (12) are 
experimental values and only ay was taken from indi- 
rect measurements. Thus it is highly expected to operate 
the thermal expansion measurement to yield a direct ay- 
The present results indicate that the range from 7^=1.83 
to 2.3 should cover the reasonable choices for the lattice 
Griineisen parameters. 

To verify these results, and also to study the behavior 
of T c as a function of pressure, we have performed explicit 
calculation based on Eq. (11). The theoretical results in 
the pressure range from to 1.0 GPa are shown p. Fig. 
2. The experimental data points of Tomita et alrB and 
Lorenz et aZJ13 measured under hydrostatic pressure con- 
ditions are also plotted for comparison. It is clearly seen 
that our calculations agree well with the experiments. 

In Fig. 3 we presented the calculated results as well as 
the experimental data points of Monteverde et aZJla and 
Deemyad et alE3 measured in the relatively high pressure 
region. Here we assume that phase transitions of all kinds 
do not occur under pressure range that we consider. We 
noticed that the experimental data points of Bccmyad et 
alE3 and the sample 4 of Monteverde et alJtB are situ- 
ated well between the two theoretical curves calculated 
by using 7g=L83 and 2.3, respectively. Interestingly, the 
agreement between our theoretical curve calculated by 
using 7g=1.83 and the experimental data points of other 
samples of Monteverde et al$3 is seen to be reasonable, 
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FIG. 3: Pressure dependence of the transition temperature 
in MgP)2 up to 40 GP&ri Experimental data am from the works 
of Monteverde et alx3 and Deemyad et alEj The inset is a 
calculation of T c under pressure up to 100 GPa. 

although there are some scatters among different sam- 
ples and the reason is not clear. Furthermore, although 
the pressure measurements are limited to the region be- 
low 33 GPa, it is seen from the inset of Fig. 3 that Eq. 
(11) continues to describe the pressure dependence of T c 
as high as 100 GPa. Even at this point, the supercon- 
ductivity is not destroyed by pressure in newly discov- 
ered superconductor MgB2. There was a discrepa*Mw-,Gtf3| 
whether pressure can destroy superconductivity.c30l23'Lll 
However, ourr-fes-ults support the conclusion of Olsen and 
collaboratorstZltJ that the possibility of destruction of 
superconductivity by the application of sufficiently high 
pressure most likely does not exist. It follows from the 
comparison of Figs. 2 and 3 that the the pressure effect 
on T c indeed provides a support to the elect ron-phonon 
mediated superconductivity in MgB2. 

In Fig. 4 we presented the normalized A and ji* as a 
function of pressure up to 30 GPa, calculated from Eq. 
(10) by using 7g=2.3. The Coulomb pseudopotential 
fi* increases slightly with pressure. Whereas A changes 
significantly with pressure. The contribution from /U*(P) 
to the variation of T c with pressure is much less important 
than that of A(P). Thus in the range from to 30 GPa 
the pressure effect of T c for MgB2 is dominated by the 
competition of A and Qd (or < uj 2 > 1 ' 2 ). 



TABLE I: Pressure dependences of superconducting state pa- 
rameters in MgB 2 . The units of d\nX/dP (X=T C , A , are 
ff c (0)) are in 10" 2 GPa" 1 . 
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FIG. 4: Pressure dependence of normalized Coulomb (fi*) 
and electron-phonon coupling (A) coupling strengths in MgB2 
calculated by using 7g=2.3. 



FIG. 5: Relation between the logarithmic volume derivative 
of N(Ef)v and the deviation £ from the full isotope effect 
exponent a = 0.5(1 — f) in nine simple sp metals supercon- 
ductors and MgB2 . The two squares are the values for MgB2 . 



Table I contains the calculated values of pressure de- 
pendences of superconducting parameters for MgB2 from 
Eqs. (5), (6), (8) and (16) by using 7G = 1.83 and 
7g = 2.3, respectively. The reliable values of j e is read- 
ily determined using Eqs. (16), (19), and (22). Wo 
obtained 7 e =1.64 and 2.07 for MgB 2 , which are close 
=1.7 for Pb and for Sn, 7 e =2.0 for Al (Ref. 
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tO 7 e 

The negative sign for dlnH c (0)/dP predicted for MgB2 
is in agreement with the measurements for all simple ele=. 
ments superconductors with the exception of thallium.Lij 
For the simple sp metals superconductors, RohrerEa has 
demonstrated that din N(Ep)v/dlnV must have act. 
proximately a value of 2.0. However, it was realizcdEa 
that the transition metals fail to show such simple be- 
havior. Our estimated din N (E F )v / dlnV =1 . 71 and 2.34 
for MgB2 are comparable to-Lhose obtained for simple 
sp metals superconduct^SEa Early measurements for 
most simple metalsa : ciaO show that there is difference 
between the quantities in dlnAo/dP and dlnT c /dP. 
This can be understood, with the aid of the results by 
Geilikman and Krcsin,E3 that is, 2A /k B T c = 3.52[1 + 
5.3{T c /uj ph ) 2 lnu ph /T c ]. The calculated data of MgB 2 
listed in Table I make it possible to support this theory. 
Since the phonon spectrum shifts under pressure, it fol- 
lows that for all superconductors with 2Ao/kBT c > 3.52 
a change of 2Ao/kBT c under pressure can be expected. It 
is interesting from the viewpoint of experiment to investi- 
gate the tunnel characteristics of MgB2 under hydrostatic 
pressure. 

Figure 5 is a plot of din N(Ep)v/dln V versus the de- 
viation £ from the full isotope effect for nine simple sp 
metals superconductors as well as MgB2. The experi- 
mental values of dlnN(Ep)v/dln V for simple metals ass 
chosen from the work of Olsen, Andres, and Geballe.El 
The experimental results for isotope effect exponent a 



and its deviation £ are taken from the works in Refs. [ 
H'fn 72j7^,^J,^5|. There are no experimental data avail- 
able for Al and In now, we_tpok the calculated results 
from Leavens and Carbotte.EJ We summarized these re- 
sults for simple sp metals in Table II. Olsen et aln sug- 
gested that din N (E p)v / dlnV is related to the isotope 
effect exponent a in metals superconductors. As seen 
from Fig. 5 the relation between din N (E p)v / din V and 
£ is not very clear when more data included. An isotope 
effect a=0.32(l) (Ref. |) in MgB 2 is consistent with what 
appears to be a systematic variation of a across the non- 
transition elements. It is well known that deviations of 
the isotope effect exponent from 1 /2 are a measure of the 
relative strengths of the Coulomb and phonon-mediated 
electron-electron interactions. It is indicated, from the 
compared values of £ for MgB2 (Ref. ^|) with Zn (Ref. 



TABLE II: Experimental values of din N(E F )v/d In V, iso- 
tope effect exponent a and its derivative £ = 1 — 2a in nine 
simple sp metals superconductors. 



Element 


Z 


din N(E F )v 

rflT, V 


a 




Zn 


2 


2.0 


0.37 


0.26 


Cd 


2 


2.9 


0.32 


0.36 


Hg(a) 


2 


1.7 


0.50 





Al 


3 


3.4 


0.325 


0.35 


Ga 


3 


1.8 


0.41 


0.18 


In 


3 


2.3 


0.466 


0.068 


Tl 


3 





0.49 


0.02 


Sn 


4 


2.3 


0.47 


0.06 


PI) 


4 


2.1 


0.478 


0.044 
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fn|, Cd (Ref. g), and Al (Ref. |§), that MgB 2 should be 
a medium coupling superconductor. 

IV. SUMMARY AND CONCLUSIONS 

The major conclusions given by present investigation 
can be summarized as follows: 

(i) A simple expression was derived for the pressure 
dependences of superconducting properties in simple sp 
superconductor on the basis of McMillan equation. The 
logarithmic volume derivatives of A, and 9d can be 
sclf-consistently determined from experiments and theo- 
ries. We gave an expression for ip from the theory of Gas-, 
pari and GyorffyO The theory of Morel and Andersoncil 
was used to obtain 0, which makes it possible to inves- 
tigate the pressure dependence of fi*. Neglecting the 
pressure dependence of n* , the present theoretical model 
can be reduced to the two popular models of Seidenc£l 
and Baryakhtar and Makarovo when taking 7 at = 2/3 
and neglecting the direct electron-electron interaction, 
respectively. Furthermore, we obtained an explicit ex- 
pression for the change of T c as a function of pressure 
with the help of Murnaghan equation. The present model 
enables us to study the pressure behaviors of some in- 
terested superconducting parameters such as the zero 
temperature energy gap Ao, the critical field at absolute 
zero temperature H c (0), the effective interaction strength 
N(Ep)v, and the electronic specific heat coefficient 7. 

(ii) We investigated the pressure effects on supercon- 
ducting properties in the newly discovered superconduc- 
tor MgB 2 using our simple approach. It was found that 
the hydrostatic pressure derivative of T c can be repro- 
duced by using the values of jg obtained from either 



the Slater relation or ab initio calculation. The calcu- 
lated d\nN(E F )v/d\n V «2.0 in MgB 2 is close to those 
obtained in simple sp superconductors. The quantita- 
tive agreement for the variation of T c with pressure in 
the low pressure region as well as high pressure region is 
very good when comparing our theoretical results with 
experimental data measured by three groups. The pre- 
dicted values of dlnH c (0)/dP, dlnAo/eLP, and 7 e are 
also comparable to those in simple sp metals supercon- 
ductors. All these characteristic pressure behaviors allow 
us to conclude that MgB 2 should be a simple electron- 
phonon mediated sp superconductor and the mechanism 
in simple sp metals superconductors is also responsible 
for the superconductivity in MgB 2 . 
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—dT c /dP observed in the low T c = 37. A ± 0.1 samples. 
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